Effects of flow fluctuations and partial thermalization on v± 
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The second and fourth Fourier harmonic of the azimuthal distribution of particles, V2 and V4,, have 
been measured in Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). The harmonic V4 
is mainly induced from V2 as a higher-order effect. However, the ratio V4/(v2) 2 is significantly larger 
than predicted by hydrodynamics. Effects of partial thermalization are estimated on the basis of a 
transport calculation, and are shown to increase «4/(«2) 2 by a small amount. We argue that the 
large value of V4/(v2) 2 seen experimentally is mostly due to elliptic flow fluctations. However, the 
standard model of eccentricity fluctuations is unable to explain the large magnitude of «4/(«2) 2 in 
central collisions. 

PACS numbers: 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

The azimuthal distribution of particles emitted in ul- 
trarelativistic nucleus-nucleus collisions at RHIC is a sen- 
sitive tool in understanding the bulk properties of the 
matter produced in these collisions (see [1] for a recent 
review). It is generally written as a Fourier series 



dN 

—— oc 1 + 2v 2 cos 2<p + 2t>4 cos ' 



(1) 



where </> is the azimuthal angle with respect to the di- 
rection of flow. In this paper, we consider analyses done 
near the center-of-mass rapidity, so that odd harmonics 
vanish by symmetry. The large magnitude of elliptic flow, 
v 2 , suggests that the lump of matter formed in a Au-Au 
collision at RHIC is close to local thermal equilibrium 
and expands as a relativistic fluid. Elliptic flow is large 
at high p t (up to 0.25 for baryons), which motivated the 
idea to study the higher-order harmonic v 4 [1, Q . Several 
analyses of V4 have been reported 0, H, @, 0] • Experimen- 
tal results give V4 ~ (v 2 ) 2 , while the ideal- fluid picture 
generally predicts V4 = \(v 2 ) 2 p|. This discrepancy has 
not yet been explained. In this paper, we investigate the 
sensitivity of V4 to two effects: viscous deviations from 
the ideal- fluid picture (Sec. IIII|). and elliptic flow fluctu- 
ations (Sec. rV]). 



II. IDEAL HYDRODYNAMICS 

We first briefly recall the prediction of relativistic hy- 
drodynamics. In this theory, the <p dependence of particle 
distribution results from a similar <p dependence of the 
fluid 4- velocity 



u(4>) = U (1 + 2V 2 cos 2(j) + 2V 4 cos ' 



■), (2) 



where <p is the azimuthal angle of the fluid velocity with 
respect to the minor axis of the participant ellipse (Io| 
(see Fig. [TJ). This is due to the fact that the overlap area 
between the two colliding nuclei is elliptic, which results 




FIG. 1: Schematic picture of a nucleus- nucleus collision de- 
picted in the transverse plane (from [0]). The principal axes 
(x' and y') of the area formed by the participants are tilted 
with respect to the reaction plane given by the axes (x and 
y) of the transverse plane. 



in anisotropic pressure gradients. For a semi-central Au- 
Au collision at RHIC, V 2 ~ 4%, and one expects V4 to 
be of much smaller magnitude, typically V4 ~ (V2) 2 . 

The fluid expands, becomes dilute and eventually 
transforms into particles. As argued in Ref. Q, fast par- 
ticles are produced where the fluid velocity is maximum, 
and parallel to the particle momentum. The resulting 
momentum distribution is a boosted thermal distribu- 
tion. Neglecting quantum statistics (this is justified in 
the transverse momentum range where V4 is measured), 
the momentum distribution for a given particle of mass 
m is 
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-p-u/T 



exp 
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(3) 



where m t = ^/p\ +m 2 , Uo(4>) = \f\ + u(cf)) 2 , and <j) is 
the azimuthal angle of the particle. Inserting Eq. @ 
into Eq. ([3]), expanding to leading order in V%, V4 and 
identifying with Eq. |T]), one obtains [8| 
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where v = U/yl + U 2 . The higher harmonic V4 is the 
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FIG. 2: (Color online) va/(v2) 2 versus pt in Boltzmann trans- 
port theory and ideal hydrodynamics for massless particles. 
Solid lines are 2-parameter fits using Eq. ([5]) over the interval 
[0.5, 2.5] GeV/c. The curves are labeled by the value of the 
Knudsen number K. Error bars are statistical. The square 
dots are results for charged pions from PHENIX averaged 
over the centrality interval 20-60%. 



sum of two contributions: an "intrinsic" V4 proportional 
to the cos 40 term in the fluid velocity distribution, V4, 
and a contribution induced by elliptic flow itself, which 
turns out to be exactly \{v2) 2 ■ The latter contribution 
becomes dominant as p t increases. 

In order to confirm these qualitative results, we solve 
numerically the equations of ideal relativistic hydrody- 
namics. The fluid is initially at rest. We choose a 
gaussian initial entropy density profile, with rms widths 
<i x = 2 fm and a y — 3 fin. The equation of state is 
that of an two-dimensional ideal gas of massless parti- 
cles, s cx T 2 , for reasons to be explained below. The nor- 
malization has been fixed in such a way that the average 
transverse momentum per particle is (pt) = 0.42 GeV/c, 
which is roughly the value for pions in a central Au- 
Au collision at RHIC [Hj]. Fig. [5] displays the varia- 
tion of Vi/(v2) 2 with the particle transverse momentum 
p t . For massless particles, m t = pt and Eq. gives 
V/t/(v2) 2 = 0.5 + k/pt, where k is independent of pt- To 
check the validity of this formula, our numerical results 
are fitted over the interval 0.5 < pt < 2.5 GeV/c by the 
simple formula 



Vi(Pt) 
V2{pt) 2 



= A + B 



<Pt) 
Pt 



(5) 



where we have introduced the average transverse momen- 
tum (p t ) in such a way that the coefficient B is dimen- 
sionless. We refer to A (resp. B) as to the induced 
(resp. intrinsic) v 4 . We find A = 0.557 and B = 0.479. 
The value of A is close to the expected value 0.5. The 
small discrepancy is due to the fact that Eqs. (|4]) are 
only valid for small values of v 2 and V4 . This approxima- 
tion breaks down at the upper end of our fitting interval, 
where «2(2.5 GeV/c) = 0.51. This large value is due 



to the fact that the equation of state is that of an ideal 
gas. For large pt, however, the intrinsic V4 term in Eq. ^ 
can be neglected, because it is linear in pt while the other 
term is quadratic in p t . Neglecting this term, the Fourier 
expansion in Eq. fl} can be done exactly. This yields 



V2n(pt) 



In(x) 



(6) 



where x — 2V2U(pt — m t v)/T, and I n (x) is the modified 
Bessel function. Inverting Eq. (|6|) with n = 1 and V2 = 
0.51, one obtains x = 1.19. Eq. ^ with n = 2 then gives 
Vi/{v2) 2 = 0.552, in better agreement with our numerical 
result. 

We have systematically investigated the sensitivity of 
our hydrodynamical results to initial conditions. With a 
smaller initial eccentricity [a x = 2 fm and a y = 2.5 fm), 
the value of A is closer to 0.5, as expected from the discus- 
sion above. We have also repeated the calculation with 
a more realistic density profile corresponding to a Au- 
Au collision at RHIC, obtained using an optical Glauber 
model calculation. We expected that B, which we un- 
derstand as the "intrinsic" V4, would be sensitive to the 
change in initial conditions, but the changes in both A 
and B were insignificant. 

Experimental results are also shown in Fig. [51 The 
value of V4/V 2 is constant, even at relatively low p t : a 
fit to these results using Eq. (O gives B = 0.01 ± 0.04, 
compatible with zero. 1 The other fit parameter is A = 
0.89 ± 0.02, significantly larger than the value 0.5 pre- 
dicted by hydrodynamics. Some of the discrepancies be- 
tween our model calculation and data can be attributed 
to the equation of state, which is much softer in QCD 
near the transition region than in our hydrodynamical 
calculation. More specifically, the coefficient B repre- 
senting the intrinsic V4 may depend on the equation of 
state. It would be interesting to investigate whether the 
small value of B seen experimentally can be attributed 
to the softness of the equation of state. On the other 
hand, our argument leading to A = i is quite general, so 
that the discrepancy with data cannot be attributed to 
the equation of state. In this paper, we investigate the 
possible origins of this discrepancy. 



III. PARTIAL THERMALIZATION 

It has been argued [l4[ that if interactions among the 
produced particles are not strong enough to produce lo- 
cal thermal equilibrium, so that the hydrodynamic de- 
scription breaks down, the resulting value of V4/(v 2 ) 2 
is higher. This is confirmed by transport calculations 
within the AMPT model [lj|. This naturally raises the 



1 Note, however, that STAR results for charged particles |13I1 
clearly display an intrinsic V4 component, although smaller than 
in our calculation. 
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FIG. 3: (Color online) Variation of the dimensionless fit pa- 
rameters A and B from Eq. © with the Knudsen number K. 
Error bars are statistical. Lines are linear fits. The points 
at K — are obtained from an independent hydrodynamical 
calculation and are excluded from the fit. 



question of how V4 reaches the hydrodynamic limit [16| . 
We investigate this issue systematically by solving nu- 
merically a relativistic Boltzmann equation, where the 
mean free path A of the particles can be tuned by vary- 
ing the elastic scattering cross section a. The degree of 
thermalization is characterized by the Knudsen number 



K = 



A 
R' 



(7) 



where R is a measure of the system size. We consider 
massless particles moving in the transverse plane (no 
longitudinal motion) [13]. In the limit K — ► 0, this 
Boltzmann equation is expected to be equivalent to ideal 
hydrodynamics, with the equation of state of a two- 
dimensional ideal gas. For sake of consistency with our 
hydrodynamical calculation, the initial phase space dis- 
tribution of particles is locally thermal: dN / d 2 xd 2 pt oc 
exp(— pt/T(x, y)), where the temperature profile T(x,y) 
is the same as in the hydrodynamical calculation. The 
Knudsen number is normalized as in Ref. [l7j : 



K = 



No 



(8) 



where N is the total number of particles in the Monte- 
Carlo simulation, and a the scattering cross section, 
which has the dimension of a length in two dimensions. 
Fig. O displays our results for two values of K. The re- 
sults for K = 0.05 are almost identical to the results from 
ideal hydrodynamics, as expected. For K = 0.5, V4,/ (v 2 ) 2 
is larger, as anticipated in Ref. [3]. Although the fit for- 
mula ([5]) is inspired by hydrodynamics, the quality of the 
fit is equally good for the Boltzmann calculation. In par- 
ticular, the ratio V4/(v2) 2 quickly saturates with increas- 
ing p t , which means that the scaling V4 oc (v 2 ) 2 still holds 
if the system does not reach local thermal equilibrium, as 
already observed in previous transport calculations [l8| . 

The sensitivity of U4 to the Knudsen number K is seen 
more clearly in Fig. [31 which displays the variation of the 
fit parameters A and B with K. A linear extrapolation 
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FIG. 4: (Color online) Results from STAR and 
PHENIX [22J for charged hadrons produced in Au-Au col- 
lisions at 200 GeV per nucleon pair, versus the number of 
participant nucleons. We have averaged the ratios Uj/(i>2) 2 
over the intervals 1.0 < pt < 2.7 GeV/c for STAR, 1.0 < 
p t < 2.4 GeV/c for PHENIX. Dash-dotted line: prediction 
from ideal hydrodynamics without flow fluctuations. Stars: 
with fluctuations inferred from the difference between t>2{2} 
and 1>2{LYZ}, Eq. (|14[) . Dotted line: eccentricity fluctuations 
from a Monte-Carlo Glauber, Eq. (|15|l . Full line: same, with 
partial thermalization taken into account, Eq. (|21[) . 



of our Boltzmann results to the limit K — gives A — 
0.524 ± 0.008 and B = 0.508 ± 0.012, to be compared 
with our results from ideal hydrodynamics A — 0.557 
and B = 0.479, in good agreement 2 . 

These transport results may be sensitive to the choice 
of initial conditions. We have assumed a locally thermal 
momentum distribution. Now, the prediction V4/(v 2 ) 2 
from hydrodynamics originates precisely from the as- 
sumption that momentum distributions are thermal in 
the rest frame of the fluid, see Eq. ([3]) . Replacing the ex- 
ponential in this equation with a more general function 
f(p ■ u) leads to v A /{v 2 ) 2 = ff"/(2f 2 ). With a Levy 
distribution f(x) — (l+x/n/T)~ n , the value of vn/{v 2 ) 2 
is enhanced by a factor (1 + n)/n. Values of n inferred 
from p t spectra of particles produced in p-p collisions are 
close to 10 [19(, which yields a slight increase from the 
prediction of hydrodynamics. 

Realistic values of the Knudsen number K, inferred 
from the centrality dependence of v 2 [20| , are in the range 
0.3 — 0.5 for semi-central collisions. For these values, 
Fig. [3] shows that U4/(«2) 2 is at most 0.6, still signifi- 
cantly below the experimental value 0.9. We conclude 
that partial thermalization alone cannot explain experi- 
mental data. 



2 There is a small residual discrepancy of a few percent between 
Boltzmann and ideal hydrodynamics, which we do not under- 
stand. 
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IV. CENTRALITY DEPENDENCE OF VM 



FLOW FLUCTUATIONS 



RHIC experiments have analyzed in detail the cen- 
trality dependence of ^/(^g) 2 - Preliminary results from 
STAR [21] and PHENIX [2^ are presented in Fig. H 
The values of Vi/(v2) 2 are larger than 0.8 for all cen- 
tralizes, and increase up to 1.6 for central collisions. 
Both experiments observe a similar centrality dependence 
of V4/(v2) 2 - STAR obtains values slightly higher than 
PHENIX. This difference may be due to nonflow effects, 
which are smaller for PHENIX than for STAR because 
the reaction plane detector is in a different rapidity win- 
dow than the central arm detector Q. Nonflow effects 
contribute both to V2 and V4. We now estimate the order 
of magnitude of the error on V4. We consider for sim- 
plicity the case when V4 is analyzed from three-particle 
correlations. The corresponding estimate of U4, denoted 
by W4{3} [23j], is defined by 



M3} = 



(008(40! - 20 2 - 20 3 )) 

M 2 



(9) 



where <pj are azimuthal angles of outgoing particles and 
angular brackets denote an average over triplets of par- 
ticles belonging to the same event. In Eq. ©, i>2 must 
be obtained from another analysis. Nonflow effects arise 
when particles 1 and 2 come from the same source [J. 
Assuming that the source flows with the same «2 as the 
daughter particles, we obtain 



(008(40! - 20 2 - 20 3 )) = v 4 (v 2 ) 2 + S nf {v 2 ) 2 



(10) 



where S n f is the nonflow correlation. The latter can be 
estimated [24| using the azimuthal correlation S pp mea- 
sured in proton-proton collisions [251 ] and scaling it down 
by the number of participants: S n f = 2S PP /N ] 



part ■ 



Di- 



viding by (V2) 4 , we obtain the corresponding error on 
Vi/{V2) 2 : 



!> 4 



2<5 D 



> 2 )V nf AWtM 2 " (n) 

In practice, the analysis is done using the event-plane 
method rather than three-particle correlations, but this 
changes little the magnitude of nonflow effects 24]. The 
error (jTTJ) varies with centrality like l/x 2 > where x ~ 
V2vN is the resolution parameter entering the flow anal- 
ysis. The numerical value 8 pp = 0.0145 has been used in 
Ref. [24[ to subtract nonflow effects from v%. It was ob- 
tained by integrating the azimuthal correlation in proton- 
proton collisions over p t . The error bar on STAR results 
in Fig. [4] is obtained using Eq. (fTTj) with S pp = 0.0145. 
The agreement with PHENIX is much improved. How- 
ever, this may be a coincidence: in the case of U4, which is 
measured at relatively large pt , nonflow effects are likely 
to be larger; on the other hand, nonflow contributions to 
V2 tend to increase V2 and decrease the ratio v^/{v2) 2 , 
which goes in the opposite direction. Finally, we must 
keep in mind that even with a rapidity gap as in the 
PHENIX analysis, there may be a residual nonflow error 
of a similar magnitude. 



The scaling V4 = 0.5 {v 2 ) 2 predicted by ideal hydro- 
dynamics only holds for identified particles at a given 
transverse momentum p t and rapidity y, for a given ini- 
tial geometry. In order to increase the statistics, how- 
ever, experimental results for v 2 and V4 are averaged 
over some of these quantities before computing the ra- 
tio V4/(v 2 ) 2 - The averaging process increases the ratio. 
For instance, the results shown in Fig. [2] are averaged 
over a large centrality interval 20-60%. Even within a 
narrow centrality class, the initial geometry varies sig- 
nificantly due to fluctuations in the initial state [26|, |27| 
We now discuss the influence of these fluctuations on v 2 
and V4. We assume for simplicity that V2 and V4 are an- 
alyzed using two-particle correlations and three-particle 
correlations, respectively. The case where the analysis is 
done using the event-plane method is more complex and 
will be discussed in Sec. IVII The estimate of V2 from 
two-particle correlations is denoted by V2{2} and defined 
by v 2 {2} 2 = (cos(20i — 202)). If v% fluctuates within 
the sample of events, (cos(20i — 202)) = ({v 2 ) 2 )- Sim- 
ilarly, if V4 and v% fluctuate, (cos(40i — 20 2 — 203)) = 
(i>4(u2) 2 )- We thus obtain 



M3} (V4(V2) 2 ) 1 {{V2Y) 



V2{2} 2 ((v 2 ) 2 ) 2 2 {{v 2 ) 2 Y 



(12) 



where, in the last equality, we have assumed that the 
prediction of hydrodynamics V4 — {V2) 2 /2 holds for a 
given value of v 2 - If v 2 fluctuates, ({v 2 ) 4 ) > ((V2) 2 ) 2 , 
which shows that elliptic flow fluctuations increase the 
observed v±/{v2) 2 - We now estimate quantitatively the 
magnitude of these fluctuations. 



A. Flow fluctuations from V2 analyses 

The magnitude of v 2 fluctuations can be inferred from 
the difference between estimates of V2, which is domi- 
nated by flow fluctuations except for very peripheral col- 
24j . The estimate from 2-particle correlations, 



lisions 



t>2{2}, gives directly ((i>2) 2 ), while the estimate of v 2 
from 4- par ticle cumulants, denoted by W2{4}, involves 

<m 4 ) m 



V2{AY^2({v 2 ) 2 ) 2 ~({v2f) 



(13) 



Inverting this relation and inserting into Eq. (|12[) . one 
obtains an estimate of the effect of V2 fluctuations on V4: 



M3} _ 1 ( M^y 

v 2 {2Y 2 I \v 2 {2}, 



(14) 



We use v{2} from [29j]; instead of «2{4}, we use the 
more recent measurement «2{LYZ} using Lee- Yang ze- 
roes [3(3, HH , which is expected to have a similar sensi- 
tivity to flow fluctuations. Data on U2{LYZ} are only 
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available for semi-central collisions. The resulting pre- 
diction for v±/(v2) 2 is shown in Fig. 2J The agreement 
with data is much improved when fluctuations are taken 
into account. We have checked numerically that our re- 
sults do not change significantly if nonflow effects are 
subtracted from ^2{2} using the parametrization intro- 
duced in Ref. H. 



B. Flow fluctuations from eccentricity fluctuations 

Since there are no data on i>2{LYZ} for the most cen- 
tral and peripheral bins, we need a model of v 2 fluctu- 
ations to cover the whole centrality range. We use the 
standard model of eccentricity fluctuations [l(| [HJ . The 
idea is that the overlap area between the colliding nuclei 
(see Fig. [IJ is not smooth: positions of nuclcons within 
the nucleus fluctuate from one event to another, even 
for a fixed impact parameter. Therefore, the participant 
eccentricity, epp, which is the eccentricity of the ellipse 
defined by the positions of participant nucleons, also fluc- 
tuates. Assuming that v 2 in a given event scales like epp, 
Eq. (H21) gives 



« 4 {3} 
V2{2} 2 



1 



=pp/ 



(15) 



c pp/ 



We estimate this quantity using the Monte-Carlo 
Glauber model [32j provided by the PHOBOS collabo- 
ration 33]. In each event, the participant eccentricity is 
defined by 



epp 



ft 



4(T 2 



(16) 



where a 2 — {x 2 } — {x} 2 and a xy = {xy} — {x} {y}, and 
{•••} denotes event- by- event averages over participant 
nucleons. Each participant is given a weight proportional 
to the number of particles it creates: 



(1 - x) + xN con / part , 



(17) 



where N co n/ palt is the number of binary collisions of the 
nucleon. The sum of weights scales like the multiplicity: 



dNch 
ch] 



pi' 



(1-x) 



part 



• xN, 



coll 



(18) 



where -/Vp art and iV co n are respectively the number of 
participants and of binary collisions of the considered 
event. We choose the value x = 0.13 which best de- 
scribes the charged hadron multiplicity observed experi- 
mentally [33)]. We define the centrality according to the 
multiplicity (fT5|) . We evaluate eccentricity fluctuations 
in centrality classes containing 5% of the total number of 
events. 

Our results are presented in Fig. 0J For peripheral 
and semi-central collisions, the estimates from eccentric- 
ity fluctuations are in good agreement with the earlier 



estimate from the difference between v 2 analyses, in line 
with the observation that this difference is mostly due to 
eccentricity fluctuations [24j . For the most central bin, 
however, eccentricity fluctuations only increase V4 L /(y 2 ) 2 
by a factor 2, while a factor 3 would be needed to match 
STAR and PHENIX data. This factor 2 can be simply 
understood. For central collisions, eccentricity fluctua- 
tions are well described by a two-dimensional gaussian 
distribution 13411: 



dN 
de 



epp 



■ exp 



pp 



c pp 
2a 2 



(19) 



This implies (e PP )/(ep P ) 2 = 2. 

We now combine the effects of flow fluctuations and 
partial thermalization, discussed in Sec. IIIII We take 
partial thermalization into account using the linear fit to 
the coefficient A from Fig. [3] 



{V2? 



2+0.18 K. 



This modifies Eq. (fT5|) into the following equation: 



M3} 
v 2 {2} 2 



+ 0.18 K 



(4p> 



-pp/ 



(20) 



(21) 



The value of K can be evaluated using the centrality 
dependence of elliptic flow. We borrow our estimates 
from Ref. 12011 . This study has recently been corrected 
and refined [2a], but the resulting estimates of K differ 
little from the original ones. Results are shown in Fig. [4] 
Partial thermalization is a small effect. Agreement with 
data is significantly improved for semicentral collisions, 
not for central collisions. For peripheral collisions, our 
calculation overshoots PHENIX data. Note that Eq. ([20]) 
was derived using the results of a Boltzmann transport 
calculation, which only applies to a dilute gas. With a 
realistic, soft equation of state, the coefficient in front of 
K could be different. 



C. A toy model of Gaussian flow fluctuations 

In order to illustrate the sensitivity of V4 to the statis- 
tics of v 2 fluctuations, we finally consider a toy model 
where the distribution of v 2 at fixed impact parameter b 
is Gaussian: 



dN 
dv 2 



1 



<7 V \/2'K 



exp 



(v 2 - Ke RP (6)) ; 
2a?, 



(22) 



where crp is the reaction-plane eccentricity obtained us- 
ing an optical Glauber model (smooth initial density pro- 
file), and k a proportionality constant. We assume that 
a v scales like N~^/ t 2 [3(|, as generally expected for ini- 
tial state fluctuations, and we adjust the proportional- 
ity constant so as to match the difference between f2{2} 
and «2{4} for midcentral collisions. The result is dis- 
played in Fig. [51 We have checked that a similar result is 
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FIG. 5: (Color online) Results using a toy model of gaussian 
v 2 fluctuations. STAR and PHENIX data as in Fig.g] Dashed 
line: ideal hydrodynamics+gaussian flow fluctuations. Full 
line: gaussian flow fluctuations and partial thermalization. 



obtained if we use the eccentricity from the Color-Glass 
condensate [13] instead of the Glauber eccentricity. For 
semicentral and peripheral collisions, this model is rea- 
sonably close to the standard model of eccentricity fluc- 
tuations. For central collisions, however, results are very 
different, because one-dimensional gaussian fluctuations 
satisfy ((U2) 4 ) / (("V2) 2 ) = 3 for central collisions, in- 
stead of 2 for eccentricity fluctuations, which are two- 
dimensional. The toy model is in very good agreement 
with data once partial thermalization is taken into ac- 
count using Eq. (T2"Tj) . However, it lacks theoretical foun- 
dations: we do not know any microscopic picture that 
would produce such gaussian fluctuations. 



FIG. 6: Effect of fluctuations on w 4 {EP}/« 2 {EP} 2 . The pa- 



rameter q, defined in Eq. (|24[1 . is plotted versus the resolution 
of the event-plane for elliptic flow. The solid curve is the usual 
case where the event-plane consists of two subevents; the dot- 
ted curve is the case where the event-plane consists of only 
one subevent [TT| ] . 



in Appendix [X] using the same methods as in Ref. [24| . 
Fig. [5] displays the variation of a with the event-plane 
resolution for elliptic flow. One sees that a < 4, which 
means that the effect of fluctuations is always smaller for 
t>4{EP} than for t>4{3}; this is confirmed by the experi- 
mental observation «4{3} > i>4{EP} [4j. The resolution 
is I when the reaction plane is reconstructed exactly. In 
this limit, w 2 {EP} = (v 2 ), v 4 {EP} = ±((v 2 ) 2 ), which 
implies a = 1. In practice, however, the maximum reso- 
lution for mid-central collisions is 0.84 for STAR [29| and 
0.74 for PHENIX 0). In the case of PHENIX, a is larger 
than 3.2 for all centralities, which means that the effect 
of fluctuations is decreased at most by 20% compared to 
our estimates in the previous section. 



VI. FLUCTUATIONS AND FLOW METHODS 



VII. DISCUSSION 



In practice, vi and V4 are analyzed using the event- 
plane method [H, [39[ . The corresponding estimates are 
denoted by W2{EP} and U4{EP}. In this Section, we ar- 
gue that flow fluctuations have almost the same effect on 
i>4{EP} as on Ui{3}. We limit our study to small fluctua- 



tions for simplicity, in the same spirit as in Ref. [24j . We 
write V2 = (v2)+5v, with (Sv) — and (Sv 2 ) = a 2 , where 
<r v characterizes the magnitude of flow fluctuations. Ex- 
panding Eq. (|I2p to leading order in a v , we obtain 



M3} = 1 , 



M2V 



«2 



Similarly, one can write 

MEP} 
v 2 {EP} 2 



1 



«2 



(23) 



(24) 



where a depends on the reaction plane resolution. A 
similar parametrization has been introduced for the fluc- 
tuations of «2{EP} II]. The expression of a is derived 



We have shown that experimental data on V4 are rather 
well explained by combining the prediction V4 = \{v2) 2 
from ideal hydrodynamics with elliptic flow fluctuations. 
If this scenario is correct, then ui/^) 2 should be inde- 
pendent of particle species and rapidity for fixed pt and 
centrality. This is confirmed by preliminary results from 
PHENIX, which give the same value for pions, kaons and 
protons [TJ]. Ideal hydrodynamics, which fails to describe 
wa(pt) for pt > 1.5 GeV/c, seems to describe well V4/V 2 
at least up to p t ~ 3 GeV/c. 

Note that our scenario does not support the picture 
of hadron formation through quark coalescence at large 
Pt [ill- We find values of V4/V2 below 1 as a result of 
the hydrodynamic expansion, which is believed to take 
place in the quark phase. But coalescence requires that 
V4,/(v2) 2 for the underlying quark distribution is much 
higher, around 2 [4l| . 

The centrality dependence of V4 offers a sensitive probe 
of the mechanism underlying flow fluctuations. Eccen- 
tricity fluctuations have been shown to explain quantita- 
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tively V2 data in Au-Au and Cu-Cu collisions. We find 
that they also explain most of the results on V4 for pe- 
ripheral and semi-central collisions. However, they are 
unable to explain the steep rise of t>4/(i> 2 ) 2 for the most 
central bins, which is clearly seen both by STAR and 

PHENIX. Data suggest that <(w 2 ) 4 ) / (M 2 ) 2 ^ 3 for 
the most central bin, while eccentricity fluctuations give 
2. Impact parameter fluctuations only increase v^jv\ by 
a few percent. We cannot exclude a priori that the large 
experimental value is due to large errors in the extrac- 
tion of V4: if we multiply the nonflow error estimated in 
Sec. |IV] by a factor 4, data agree with our calculation 
for central collisions; however, the agreement is spoilt for 
peripheral collisions. It therefore seems unlikely that the 
discrepancy is solely due to nonflow effects. These re- 
sults suggest that initial state fluctuations do not reduce 
to eccentricity fluctuations, as recently shown by a study 
of transverse momentum fluctuations [42| . Interestingly, 
the direct measurement of v 2 fluctuations attempted by 
PHOBOS [HI, which agrees with the prediction from ec- 
centricity fluctuations, does not extend to the most cen- 
tral bin. 

An independent confirmation that ((v 2 ) 4 ) / ((v 2 ) 2 ) — 
3 for central collisions could be obtained from the 4- 
particle cumulant analysis. Interestingly, there is no 
published value of t> 2 {4} for the most central bin: the 
reason is probably that v 2 {4} cannot be defined using 
Eq. (fT3"j) , because the right-hand side is negative. This in- 
dicates that ((V2) 4 ) I ((V2) 2 ) > 2. It would be interest- 
ing to repeat the cumulant analysis for central collisions, 
and to scale the right-hand side of Eq. ([M]) by v 2 {2} 4 . 

The ratio should be around —I if ((v 2 ) 4 ) / ((w 2 ) 2 ) 2 ^ 3. 
This would give invaluable information on the mechanism 
driving elliptic flow fluctuations. 



resolution, the last term is the contribution of fluctua- 
tions to v 2 {EP}. The definition of a, Eq. can be 
rewritten as 



S(V4/(V2) 2 

(V(« 2 ) 2 



(V2Y 



(A3) 



The three terms in Eq. (|A2|) give additive contributions 
to a, which we evaluate in turn. 

We start with the correlation with the event-plane. 
The event plane is determined from elliptic flow [38| . 
Even flow harmonics w 2n are analyzed by correlating 
particles with this event plane: (cos2n(0 — \&i 
v 2n7t2n(x)i where the resolution 7\L 2n is given by [3£ 



K 2n {x) = — e 



-x 2 /2 




(A4) 



where x is the resolution parameter, which is estimated 
using the correlation between two subevents. For n = 2, 
this equation reduces to 



1 + x 2 



x z 



(A5) 



These relations are derived neglecting flow fluctuations. 
If V2 fluctuates, the resolution parameter \ scales like w 2 , 
X = fV2- Assuming in addition that V4 scales like (f 2 ) 2 , 
the relative change due to fluctuations is, to leading order 



*<cos4fa-¥ fl )) _ a 2 v {{v2) 2 n A (rv 2 )) 



(cos4(>- 



2 (V2) 2 n 4 (r (v 2 )) 



2(«2> 



(A6) 



APPENDIX A: EFFECT OF FLUCTUATIONS ON 
THE EVENT-PLANE v 4 

In this Appendix, we derive the expression of a in 
Eq. (|24| . This parameter measures the effect of fluctua- 
tions on i>4/(v 2 ) 2 when flow is analyzed using the event- 
plane method. The event plane V4 is defined by 



v A {EP} 



(cos4(</>- 

Ra 



(Al) 



where <p is the azimuthal angle of the particle, '5 r is the 
angle of the event plane, and i?4 is the event-plane resolu- 
tion in the fourth harmonic. Using Eq. (|A1|) . the relative 
variation of V4/(i> 2 ) 2 due to eccentricity fluctuations can 
be decomposed as the sum of three contributions 

<WM 2 ) = *(cos4(^-tt fl )) _S_R4 _ 2 S_V2 (A2) 
(VM 2 ) (cos 40/)-^)) R A v 2 ' 1 ' 

The first term on the right-hand side is the contribution 
of fluctuations to the correlation with the event plane, 
the second term is the contribution of fluctuations to the 



where the right-hand side is evaluated for x = r ( v 2), 
the average resolution parameter. Using Eq. (|A5|) . one 
obtains 



f 



(X ^4 (x)) = 



2x 2 - 1) 



l + e x 2 (x 2 -l) 



(A7) 



Inserting into Eqs. (|A6[) and (|A2[) , and identifying with 
Eq. (|A3[) . we obtain the contribution to a from the cor- 
relation with the event plane: 



X 2 (e* +2 X 2 -1) 
f+ e x 2 ( x 2_ 1) ' 



(A8) 



We now evaluate the second term in Eq. (|A2[) , namely, 
the shift in the resolution from fluctuations. The reso- 
lution is defined as R4 = 7?4(x exp ), where x cxp is deter- 
mined from the correlation between subevents. Flow fluc- 
tuations shift the estimated resolution. Writing x° xp = 
X + ^X? one obtains, to leading order in 5x, 



SR 4 = xK(x) s x 

R4 7^4 (x) X ' 



(A9) 
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Eq. (|A5|) gives 



where /q.i is a shorthand notation for Jo,i(x 2 /2)- 



xK( x ) _ 2( e x - x 2 -i) 
n4x) l + ex 2 (x 2 -l)' 



(A10) 



The shift in the resolution to fluctuations is given by 
Eq. (A7) of Ref. [H 



S_X 

X 



2(v) 



1-2X 2 



(All) 



where io,i is a shorthand notation for io,i(x 2 /2), and Xs 
denotes the resolution parameter of a subevent: Xs — 
x/V% in the usual case when the event plane consists of 
two subevents (38|, and Xs = X if the event plane has 
only one subevent Inserting Eqs. (j A10[) and (|A11[) 
into (K9\i and (|A"2)) . and identifying with Eq. ([A3]) , we 
obtain the contribution to a from the resolution: 



X 



1 



l-2xf. 



4/7 



2) is 



^ l + e x 2 ( x 2_! 
Finally, the third term in Eq. 

2 = — — 2 ("«2 - !) 

f 2 (« 2 ) 2 

where a„ 2 is given by Eq. (23) of Ref. [24| 

In - h /_ o _ o 4if 



= 2- 



' '2x 2 -2 X 2 



(A12) 



(A13) 



(A14) 



The final result is obtained by summing the three con- 
tributions from Eqs. (|A8j) . (|A"I2|l and ifAli} : 



1). 



(A15) 



The limit of low resolution x — > (resp. high resolution 



X — > oo) is a ep = 6 (resp. 1) 
(resp. 1), a = 4 (resp. 1). 



1 (resp. 0), a v 
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